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We assume that due to their small binding energy, 
deuterons from the midrapidity region cannot survive in¬ 
side fireballs created in ultrarelativistic nuclear collisions. 
All observed deuterons are thus expected to be created 
from a neutron proton pair at freeze-out by coalescence. 
This coalescence process involves small relative momenta 
and is thus a non-relativistic phenomenon. Switching to 
center-of-mass and relative coordinates and momenta, the 
probability dVd(i?, Pd) for creating a deuteron with given 
momentum Pd at space-time point R is most conveniently 
evaluated in the center-of-mass frame of this deuteron. 
The highly relativistic dynamics of the hreball needs to 
be fully taken into account only later on, when calculat¬ 
ing the Lorentz invariant deuteron spectrum by integrating 
over the freeze-out hypersurface cr: 

Pd^ = /Pd • dV(P) diVd(P,Pd) (1) 

In non-relativistic statistical quantum mechanics dN^ 
is calculated by projecting the deuteron density matrix 
on the density matrix of the fireball at freeze-out Q. 
The latter contains, however, the relativistic dynamics of 
the heavy ion collision and is difficult to compute. This 
problem may be by-passed by switching to the equivalent 
Wigner function formalism: 

P. Danielewicz has shown in how, starting from a 
complete quantum mechanical description in terms of den¬ 
sity matrices, one can via a set of reasonable assumptions 
reduce the deuteron formation rate to an integral over a 
product of nucleon single particle Wigner functions (which 
we will approximate by hydrodynamical, local equilibrium 
distribution functions) and a quantum mechanical transi¬ 
tion matrix element. The transition matrix element im¬ 
plements conservation of energy-momentum in deuteron 
formation via interactions with third particles in the fire¬ 
ball and determines the probability of deuteron formation. 

We derive the following expression for the number dN^ 
of deuterons with four-momentum Pd freezing out at 
space-time point R: 

diVd(P,Pd) = 7:^ /d(P,Pd) H^R) R(P,Pd) , (2) 

[Ztt 

where the first term accounts for spin degeneracy and 
phasespace, /d is the local thermal equilibrium distribution 
function for an elementary particle with deuteron quantum 
numbers, and H describes the dependence of the fugacity 
exp{fj./T) (and thus of the particle density) on the posi¬ 
tion in the fireball. A implements the quantum mechanics 
of the formation process and is determined by the inter¬ 
nal structure of the deuteron and the variation of the nu¬ 
cleon single particle distribution functions / around the 
deuteron formation point R: 

R(P,Pd) = J V{r,q) X 


U{f=0,q^0) P2(r = 0) 

Here r and q are the relative coordinates and momenta 
of the neutron and proton in the deuteron center-of-mass 
frame; V is the Wigner transform of the internal wave- 
function ipd of the deuteron 

P(r, J cfx ip*^{f+x/2)(pd{r-x/2) . (4) 

The variation of the local equilibrium distribution func¬ 
tions / is due to the hydrodynamic expansion of the fire¬ 
ball. A is sensitive to flow gradients in very much the 
same way as the correlation function in two-particle in¬ 
terferometry iQ. In the absence of flow gradients A = 1 
for sufficiently large hreball volumes; for smaller volumes 
one recovers the model proposed by Hagedorn for pp col¬ 
lisions who suggested A = /q dV as the probability 
for deuteron production in a small reaction volume fl 0. 

Experimental values for the invariant coalescence factor 
B 2 = d/p^ cluster around (1 — 2) • 10“^ GeV^ for heavy 
ion collisions at the BEVALAC, and are up to a factor 10 
smaller in newer experiments at the AGS and SPS. 

Most recent measurements for B 2 come from the GERN 
Newmass/NA52 ||^ experiment for ^°®Pb-|-Pb collisions at 
158GeV/A: B 2 = (1.3 ± 0.2) • IQ-^GeV^ for deuterons 
and B 2 = (1.1 ± 0.3) • 10“^ GeV^ for antideuterons. To 
compare numerical results of our model with these exper¬ 
imental data, hreball size and how parameters for the hy¬ 
drodynamical description of the reaction zone are needed. 
In our group, efforts are presently being made to deter¬ 
mine them from single particle and pion correlation spec¬ 
tra. It should be noted that except for those parameters 
our model is parameter free; the deuteron spectrum fol¬ 
lows self-consistently from other spectra. First estimates 
for the NA52 data have shown to give the correct order of 
magnitude for B 2 , with values for A of ^ 0.6 — 0.75 at the 
space point with maximal thermal emission. 
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